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Line formation in solar granulation: 
V. Missing UV-opacity and the photospheric Be abundance 
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Abstract. The possibility of unaccounted for opacity sources in the UV for late-type stars has often been invoked to explain 
discrepancies between predicted and observed flux distributions and spectral line strengths. Such missing UV-opacity could 
among other things have a significant impact on abundance determination for elements whose only relevant spectral features 
are accessible in this wavelength region, such as Be. Here, the study by Balachandran & Bell (1998) is re-visited in the light 
of a realistic 3D hydrodynamical solar model atmosphere and the recently significantly downward revised solar O abundance 
obtained with the same model atmosphere. The amount of missing UV-opacity, if any, is quantified by enforcing that the OH 
A-X electronic lines around 313 nm produce the same O abundance as the other available diagnostics: OH vibration-rotation 
and pure rotation lines in the IR, the forbidden [O I] 630.0 and 636.3 nm lines and high-excitation, permitted O I lines. This 
additional opacity is then applied for the synthesis of the Be II line at 313.0nm to derive a solar photospheric Be abundance in 
excellent agreement with the meteoritic value, thus re-enforcing the conclusions of Balachandran & Bell. The about 50% extra 
opacity over accounted for opacity sources can be well explained by recent calculations by the Iron Project for photo-ionization 
of Fel. 
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1. Introduction 

Among other things, the light elements Uthium, beryllium and 
boron offer an opportunity to extract information about stel- 
lar interiors and evolution. Due to their fragile nature, these 
elements are destroyed by nuclear processing when brought 
to sufficiently high temperatures by stellar convection or other 
mixing events. Since the destruction of Li, Be and B occurs at 
somewhat different temperature regimes 2.6 • 10^ K for ^Li, 
3.5 • 10'' K for Be and 5 • 10'^ K for B), the relative amount 
of depletion of these elements can function as sensitive probes 
how deep the stellar mixing has proceeded. Indeed there is ex- 
tensive information available on light element depletion in stars 
of different spectral types (e.g. Deliyannis et al. 2000), con- 
firming the presence of additional mixing over that predicted 
by the standard mixing length theory for convection during the 
pre-main sequence and main sequence evolution. Among the 
invoked explanations for the observed depletion patterns, slow 
rotationally-induced mixing has perhaps received the most at- 
tention recently given the reasonable agreement with observa- 
tions for its predictions over a range of stellar environments 
(DeHyannis et al. 2000). 

Neither has our Sun been immune to light element deple- 
tion. The observed photospheric Li abundance is about a factor 



of 140 lower than measured in Cl-chondrite meteorites (Miiller 
et al. 1975; Kiselman 1997; Asplund et al. 2003a; Lodders 
2003), a difference which is far beyond the uncertainties in- 
volved in the photospheric and meteoritic abundance deter- 
minations. Until recently, it was believed that also Be and B 
were depleted in the solar photosphere by factors of about 1.8 
(Chiemlewski et al. 1975) and 1.9 (Kohl et al. 1977; Kiselman 
& Carlsson 1996), respectively, relative to the meteoritic evi- 
dence. This standard picture was challenged by Balachandran 
& Bell (1998, see also Bell et al. 2001), who attempted to 
calibrate the since-long suspected missing UV-opacity (e.g. 
Vernazza et al. 1976; Gustafsson & Bell 1979) by enforcing 
that the same oxygen abundance should be determined from the 
OH A-X electronic lines around the 313 nm region, where the 
crucial Be II lines are located, as the more reliable IR vibration- 
rotation counterparts of OH (log eo = 8.91' in the case of the 
ID Holweger-Miiller solar model atmosphere). Their investi- 
gation indicated a missing opacity for these wavelengths of 
about 60% over known opacity sources, which they attributed 
to photo-ionization of Fe I. Indeed, subsequent analysis of de- 
tailed atomic calculations performed within the framework of 
the Iron Project (Hummer et al. 1993) indicated that this hy- 
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pothesis was quite plausible (Bell et al. 2001). When account- 
ing for this missing continuous opacity in the synthesis of the 
Be II doublet at 313 nm, Balachandran & Bell (1998) found that 
the resulting photospheric Be abundance was indistinguish- 
able from the coiimionly adopted meteoritic Be abundance of 
logeBo = 1-41 (Lodders 2003). Furthermore, both the pho- 
tospheric and meteoritic B abundances have recently been re- 
vised, bringing them into apparent agreement within their (sub- 
stantial) uncertainties (Cunha & Smith 1999; Zhai & Shaw 
1994). Current thinking thus implies that while Li in the so- 
lar convection zone has been significantly destroyed, neither 
Be nor B have been depleted relative to their proto-solar val- 
ues. As indicated above, such a depletion pattern requires addi- 
tional non-standard mixing, although the precise nature of this 
mixing is still eluding astronomers. 

It is important to reahse that the influential analysis of 
Balachandran & Bell (1998) hinges on at least six crucial as- 
sumptions: 

1. the standard ID hydrostatic solar model atmospheres em- 
ployed in the analysis provide a realistic description of the 
line formation, 

2. the adopted transition probabiUties for the OH lines are cor- 
rect, 

3. the adopted microturbulence is appropriate for these partly 
saturated OH A-X Unes, 

4. the continuum placement in the very crowded UV region of 
the OH lines is correct, 

5. local thermodynamic equilibrium (LTE) for the OH line 
formation and equilibrium chemistry for the OH molecule 
formation are justified, and 

6. LTE is good approximation for the Ben line formation. 

Addressing some of these lingering potential sources of errors 
is now possible. Recently, the new generation of 3D hydrody- 
namical model atmospheres (Stein & Nordlund 1998; Asplund 
et al. 1999, 2000a; Asplund & Garcia Perez 2001) has been 
applied to studies of spectral line formation in the solar atmo- 
sphere for abundance determinations (Asplund 2000; Asplund 
et al. 2000b, 2003a,b,c,d; Allende Prieto et al. 2001, 2002b; 
Shchukina & Trujillo Bueno 2001). A particularly noteworthy 
finding from these studies is the substantial downward revision 
of the solar photospheric O abundance to log eo = 8.66 ± 0.05 
(Asplund et al. 2003d). This new abundance has been deter- 
mined from a multitude of different spectral lines and for the 
first time all O diagnostics yield consistent results: the forbid- 
den [Ol] line at 630 nm, high-excitation, permitted Ol, OH 
vibration-rotation fines and OH pure rotation lines. In view of 
the significantly lower O abundance commonly accepted today 
(Lodders 2003) compared with the value believed to be cor- 
rect at the time of the study by Balachandran & Bell (1998), it 
is clearly of importance to repeat their analysis using the new, 
highly reafistic, 3D hydrodynamical solar model atmosphere. 
Such calculations are presented here, basically confirming the 
findings of Balachandran & Bell (1998) of significant (~ 50%) 
missing UV-opacity and a photospheric Be abundance in agree- 
ment with the meteoritic abundance. Of the above-mentioned 
six uncertainties in the original Balachandran & Bell (1998) 



study, apparently now only the penultimate one remains which 
could possibly modify these conclusions. 

2. 3D spectral line formation calculations 

The same 3D hydrodynamical solar model atmosphere which 
has previously been confronted successfully with a range of 
observational diagnostics (e.g. Asplund et al. 2000a) and used 
for abundance analysis purposes (Asplund 2000; Asplund et 
al. 2000b, 2003a,b,c,d; Allende Prieto et al. 2001, 2002b; 
Shchukina & Trujillo Bueno 2001) is employed here. For nu- 
merical details of the simulation, the reader is referred to Stein 
& Nordlund (1998), Asplund et al. (2000a) and Asplund & 
Garcia Perez 2001). 

The 3D spectral line formation calculations are performed 
under the assumption of LTE for the ionization and exci- 
tation balances and for the source function {Si, = B^). 
Instantaneous chemical equilibrium is furthermore assumed for 
the OH molecule formation. The radiative transfer is solved 
for 17 different inclined directions (AT^ = 7V<^ = 4 plus 
the vertical) using realistic background continuous opacities 
(Gustafsson et al. 1975; Asplund et al. 1997) and an equation- 
of-state which accounts for excitation, ionization and molecule 
formation of the most important elements (Mihalas et al. 1988). 
The photo-ionization cross-sections of Fe I as calculated by the 
Iron Project (Hummer et al. 1993) have not been included. It is 
likely that the missing UV-opacity quantified below is indeed 
due to Fei bound-free opacity (Bell et al. 2001). The inclined 
rays are disk-integrated adopting a solar rotational velocity of 
1.8kms~^ to produce flux profiles; the exact choice of rota- 
tional broadening is however inconsequential for this study. 
The instantaneous flux profiles are subsequently averaged over 
a simulation time-sequence corresponding to about 50 min of 
solar time, which is sufficiently long to yield statistically sig- 
nificant results. 

A major advantage with the new generation of 3D hydrody- 
namical model atmospheres is that the traditional free param- 
eters of stellar spectroscopy (mixing length parameters to de- 
scribe convection and the micro- and macroturbulence for line 
broadening) no longer are needed due to the self-consistently 
calculated convective velocities and the corresponding Doppler 
shifts (Asplund et al. 2000a,b; Asplund 2000; Allende Prieto et 
al. 2002a). Thus, no microturbulence ^mrb enters the 3D calcu- 
lations presented herein, although the corresponding ID calcu- 
lations have been performed with both ^turb ~ l.Okms^^ (the 
standard case as well as the value adopted by Balachandran & 
Bell 1998) and ^turb = 1.2 km s^^, to assess the uncertainties 
attached to the adopted microturbulence in the original analysis 
of Balachandran & BeU (1998). 

The molecular data for the OH and Be lines are well- 
known. Here a small selection of four OH A-X lines are cho- 
sen (OH 312.80, 312.82, 313.91 and 316.71 nm) covering flie 
relevant range of excitation potential. As will be evident later 
on, all fines suggest very similar multiplication factors to the 
opacity at these wavelengths, and the exact choice of OH lines 
is thus unimportant. The adopted ^/-values for the OH A-X 
lines are taken from GilUs et al. (2001), which are based on the 
most recent theoretical calculations and laboratory measure- 
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Table 1. The predicted line strengths in pm (=10 mA) for the OH A-X and Be Hnes using ID MARCS and Holweger-Miiller 
solar model atmospheres calculated with different combinations of microturbulence (^turb) and additional continuous opacity 
over standard values (Ak™"*/kJ;°"*). The "observed" equivalent widths used for the 3D calculations are those computed using 
Ak™"7k™"* = 0.60 (i.e. 60% extra opacity) and ^turb = 1-0, as found by Balachandran & Bell (1998). 



line 




log 3/ A«:r'V«:r'' (MARCS) 


Wx (Holweger-Miiller) 


[nm] 


[eV] 


logeo =8.75 


logeo = 8.91 






logeBo = 1.40 


logeBc = 1.40 






Cturb ~ 1.0 Cturb ~ 1.2 


Cturb = 1.0 Cturb = 1.2 






[kms^M [kms-1] 


[kms-^] [kms-^] 



OH 312.80 0.102 

OH 312.82 0.209 

OH 313.91 0.760 

OH 316.71 1.108 

Ben 313.10 0.000 



2.425 0.00 
0.60 

2.074 0.00 
0.60 

1.563 0.00 
0.60 

1.544 0.00 
0.60 

0.468 0.00 
0.60 



6.53 


6.69 


5.90 


6.04 


7.55 


7.75 


6.89 


7.07 


7.47 


7.67 


6.78 


6.95 


6.15 


6.29 


5.46 


5.58 


8.79 


8.91 


7.63 


7.73 



6.55 


6.71 


6.01 


6.15 


7.58 


7.78 


7.02 


7.21 


7.51 


7.71 


6.92 


7.10 


6.19 


6.33 


5.59 


5.72 


8.63 


8.74 


7.63 


7.73 



ments. For completeness, it is noted that a dissociation energy 
of 4.392 eV (Sauval & Tatum 1984) has been employed but 
it has no influence on the conclusions obtained here. The Be II 
313.10nm5/-value comes from the VALD database (Piskunov 
et al. 1995). The relevant input data is listed in Table[2 

The study presented here is done strictly differentially be- 
tween the ID and 3D cases, thereby avoiding a direct compar- 
ison with observations. Thus, in effect an intrinsic assumption 
for these calculations is that the investigation of Balachandran 
& Bell (1998) is done correctly in terms of the comparison with 
observations within the framework of ID model atmospheres 
and line formation. The main reason to here avoid confronta- 
tion with observations is that the current version of the 3D line 
formation code can only treat lines of one element or molecule 
at a time, while the OH A-X and Be lines are located in a very 
crowded spectral region with contributions from many different 
species. An advantage with the procedure adopted here is that 
the impact of the 3D model atmospheres and the new low so- 
lar O abundance can be isolated without the possible confusion 
introduced by for example continuum placement and choice of 
transition probabilities. As there are all reasons to believe that 
the careful analysis of Balachandran & Bell (1998) was done 
correctly, we are still in a position to evaluate the solar photo- 
spheric Be abundance and the amount of missing UV-opacity. 

As the predicted line shapes are intrinsically different for 
the ID and 3D cases even for the same line strengths, the 
here adopted procedure is simplified to only use equivalent 
widths. With the absence of direct observational confrontation, 
the observed line strengths of the relevant OH A-X lines have 
been estimated by redoing the analysis of Balachandran & Bell 
(1998) with the continuous opacities multiplied by a factor of 
1.6. The adopted solar O abundances for the calculations us- 
ing the ID MARCS (Asplund et al. 1997) and the Holweger- 
Miiller (1974) solar model atmospheres are log eo = 8.75 and 
log eo = 8.91, respectively, as found by Balachandran & Bell 



(1998) from the OH IR vibration-rotation lines with the two 
model atmospheres. In addition, the same calculations have 
been performed without the additional continuous opacity as 
well as for the two values of the microturbulence. The pre- 
dicted ID equivalent widths of the OH A-X lines for the dif- 
ferent cases are listed in Table [J together with the correspond- 
ing results for the Be II 313.10nm line. The equivalent widths 
computed with a multiplication factor of 1.6 to the continuous 
opacity will thus function here as substitutes for the observed 
line strengths. 

From a closer inspection of Table [2 it is clear that the 
two ID model atmospheres do not yield exactly the same 
"observed" equivalent widths, even if Balachandran & Bell 
(1998) claimed that their analysis yielded the same multipli- 
cation factor to the continuous opacities in the two cases. 
The reason for this relatively small discrepancy has not been 
identified but can probably be traced to a combination of 
slightly different adopted input data (continuous opacities etc), 
model atmospheres (exact version of MARCS models, pressure- 
integrated Holweger-Miiller model atmosphere or not), overall 
solar chemical composition, and numerical implementation of 
the ID spectrum synthesis. Assuming for the moment that the 
calculations with the Holweger-Miiller model atmospheres are 
correct with 60% extra opacity, the same result would be ob- 
tained for the MARCS case with about 50% additional opac- 
ity instead. It should be emphasized that this difference has no 
significant impact on the conclusions regarding missing UV- 
opacity and the photospheric Be abundance. 

3. Results 

Table ^ and Table |2] list the predicted line strengths for the 
OH and Bell lines using the ID and 3D model atmospheres, 
respectively. As described above, the theoretical equivalent 
widths calculated with ID model atmospheres, a microturbu- 



4 



Asplund; Line formation in solar granulation: V. The photospheric Be abundance 



1.2 




0.2 - 



0.0 1 ' 

312.795 312.800 312.805 3 1 2.81 3 1 2.8 1 5 312.820 
Wavelength [nm] 

Fig. 1. The dotted line denotes the OH 312.80nm Une calcu- 
lated with the MARCS model atmosphere using a solar oxygen 
abundance of logeo = 8.75, a microturbulence of ^tmb ~ 
1.0 km s^^ and 60% additional continuous opacity over the 
standard value computed by the spectrum synthesis program. 
According to Balachandran & Bell (1998) this fits the observa- 
tions very well, and is here used as a substitute to the observed 
profile in this crowded spectral region. The dashed line corre- 
sponds to a 3D line calculation with log eo ~ 8.66 and without 
any extra opacity. The solid line represents the 3D results with 
the same abundance but with 50% extra continuous opacity, 
which yields the same overall line strength as the here shown 
ID profile. No micro- and macroturbulence enter the 3D calcu- 
lations but the Doppler shifts caused by the convective motions 
make the line broader than the ID case which is here shown 
without macroturbulence. 

lence of ^tuib = 1.0 km s^^ and a multiplication factor of 1.6 
to the continuous opacities around 313 nm will here function 
as the "observed" values to be reproduced with the predictions 
from the 3D spectrum synthesis. For the OH lines in 3D a solar 
O abundance of log eo = 8.66 has been adopted, as indicated 
from our recent analysis of the OH vibration-rotation, the OH 
pure rotation, the forbidden [O l] lines at 630.0 and 636.3 nm, 
and high-excitation, permitted O I lines (Asplund et al. 2003d). 
It is then straightforward to compare the predicted equivalent 
widths using the 3D hydrodynamical model atmosphere with 
different multiplication factors to the continuous opacity with 
the "observed" line strengths of the OH A-X lines. It is imme- 
diately obvious that without additional continuous opacity over 
those included in the spectrum synthesis package the OH Unes 
are significantly too strong (Fig. [Q. Using an average of the 
"observed" values from the Holweger-Miiller and MARCS ID 
model atmospheres, this missing opacity can be estimated to 
be about 50%, which is indeed very similar to the value of 60% 
obtained by Balachandran & Bell (1998). 

The "observed" equivalent width of the Be II 313.0nm line 
is independent of whether the Holweger-Miiller or the MARCS 
ID model atmosphere have been used to estimate it, as clear 
from Table n (7.63 pm=76. 3 mA). In 3D, without any addi- 
tional continuous opacity this line strength would imply a solar 
Be abundance of log ebc = 1-23. However, as clear from the 
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313.095 313.100 313.105 313.1 10 313.1 15 313.120 
Wavelength [nm] 

Fig. 2. The dotted line denotes the Bel 313.10nm line calcu- 
lated with the MARCS model atmosphere using a solar beryl- 
lium abundance of log ebc — 1-40, a microturbulence of 
fturb = l.Okms^^ and 60% additional continuous opac- 
ity, which is required to fit the observations according to 
Balachandran & Bell (1998). The dashed line corresponds to 
a 3D line calculation with log eec = 1-38 and without any ex- 
tra opacity. The solid line represents the 3D results with the 
same abundance but with 50% extra continuous opacity, which 
yields the same line strength as the here shown ID profile. 

comparison of the OH A-X lines, there are strong indications 
that the missing UV-opacity amounts to about 50%. Taking 
this extra opacity into account the 3D spectrum synthesis leads 
to accordingly weaker lines, yielding a higher solar Be abun- 
dance of log ebg = 1-38 (Fig.|5J. Combined with an estimate 
of the errors in the missing opacity and the overall agreement 
between predicted and observed line profiles (Balachandran & 
Bell 1998), the final result is thus log ebo = 1.38 ± 0.09. This 
best estimate using a 3D model atmosphere is indistinguishable 
from the meteoritic value of log eBc = 1-41 ± 0.08 (Lodders 
2003) given the uncertainties in both the photospheric and me- 
teoritic analyses. This is even more true in the light of the pro- 
posed downward revision by 0.04 dex to the meteoritic abun- 
dance scale caused by the re-analysis of the photospheric Si 
abundance on which the absolute meteoritic abundance scale 
is anchored (Asplund 2000), which would bring the meteoritic 
value down to log ebg = 1-37 ± 0.08. 

The findings of Balachandran & Bell (1998) of significant 
missing UV-opacity and no Be depletion within the solar con- 
vection zone throughout the 4.6 Gyr of solar evolution until to- 
day is therefore confirmed by this analysis based on a realistic 
3D model atmosphere rather than ID models. 

4. Remaining uncertainties 

In this section, uncertainties which may have a bearing on the 
findings are discussed. 

Model atmospheres: Available evidence from various con- 
frontations between the predictions from the here employed 
3D solar simulation with observed diagnostics all suggest that 
indeed the new 3D hydrodynamical model atmospheres are 
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Table 2. The predicted line strengths in pm for the OH A-X and 
Be Hnes using a 3D hydrodynamical solar model atmosphere. 
The OH lines are computed for log eo = 8.66 and the Be line 
forlogeBc — 1-38. 



line Xcxc iog gf W\ Wx 

[nm] [eV] A«:^°"7«:^°"' 

= 0.00 = 0.50 

OH 312.80 0.102 -2.425 6.38 5.90 

OH 312.82 0.209 -2.074 7.42 6.91 

OH 313.91 0.760 -1.563 7.29 6.75 

OH 316.71 1.108 -1.544 5.87 5.34 

Ben 313.10 0.000 -0.468 8.62 7.65 



highly realistic, representing a significant improvement over 
existing ID models. It is reassuring that the main conclusions 
regarding the missing UV-opacity and the solar Be abundance 
are robust against choice of model atmosphere, whether in ID 
or in 3D. 

Line broadening: It should be noted that the line strengths 
of the OH A-X lines are somewhat sensitive to the adopted mi- 
croturbulence in ID. Had instead Balachandran & Bell (1998) 
employed a larger microturbulence, they would have obtained 
slightly smaller multiplication factors to the continuous opac- 
ity, as clear from Tabled Due to the differential nature of this 
study relative to that of Balachandran & Bell (1998), the mul- 
tiplication factor estimated here using a 3D solar model atmo- 
sphere would then have been adjusted slightly accordingly. For 
reasonable values of the microturbulence, however, the changes 
are small to the estimated amount of missing UV-opacity and 
the derived solar Be abundance. Available solar analyses all in- 
dicate that a microturbulence of 1 . ± . 1 km s ^ ^ for flux spec- 
tra is appropriate when relying on ID model atmospheres. The 
uncertainties attached to the exact choice of microturbulence 
can therefore in this case be safely ignored. 

Transition probabilities: Conceivably, the adoption of er- 
roneous transition probabilities for either the OH vibration- 
rotation or the OH A-X lines could be misleading, yielding in- 
correct estimates of the amount of missing UV-opacity, if any, 
and consequently the derived solar Be abundance. Fortunately, 
the necessary ^/-values for both sets of lines are accurately 
known from theoretical calculations and laboratory measure- 
ments (Goldman et al. 1998; Gillis et al. 2001), which are very 
similar to those used by Balachandran & Bell (1998) in their 
study. The uncertainties in the transition probabilities are suf- 
ficiently small not to compromise the conclusions presented 
here. Similarly, as long as consistent partition function and 
molecular equilibrium constants are employed, the results are 
independent of such input data. 

Continuum placement: A potentially more serious short- 
coming of the study of Balachandran & Bell (1998) and thus 
also the present analysis, is the exact placement of the contin- 
uum level. If the observed continuum is estimated to be too 
low, a too low abundance is derived, in particular given the fact 
the relevant lines for these studies are at least partly saturated. 
The region around the OH A-X and Be lines at 3 1 3 nm is noto- 
riously crowded and it is difficult to identify wavelength win- 
dows which traces the continuum. The problem is compounded 



by the fact that many of the lines in this region remain uniden- 
tified, rendering even the use of spectrum synthesis to locate 
small continuum portions uncertain. As this study does not 
involve direct comparison with observations, it relies entirely 
on that this part in the original Balachandran & Bell (1998) 
work was carried out properly. From the relatively small wave- 
length regions shown in their figures it appears however that 
the continuum placement is essentially correct, which should 
minimize the uncertainties for the final results. 

Non-LTE line formation for OH lines: While LTE is in all 
likelihood an excellent approximation for the formation of the 
OH vibration-rotation lines, this is not guaranteed for the OH 
A-X lines (Hinkle & Lambert 1975). To date, very little work 
on non-LTE line formation has been invested for molecules 
with the notable exception of CO lines (e.g. Uitenbroek 2000). 
For OH A-X lines, the author is only aware of the very prelim- 
inary study of Asplund & Garcia Perez (2001) based on only a 
two-level OH model. A fundamental obstacle in this respect is 
the lack of data for excitation and ionization by collisions with 
electrons and hydrogen atoms, even more so than for the cor- 
responding case of atoms. With their very simplified approach 
Asplund & Garcia Perez (2001) indeed found the possibility of 
significant departures from LTE for the OH A-X lines, amount- 
ing to about 0.2 dex in terms of oxygen abundance for the so- 
lar case in ID. If this significant weakening of the predicted 
OH A-X lines in non-LTE is correct, the conclusions regard- 
ing missing UV-opacity must be revised. Indeed, a non-LTE 
effect of 0.2 dex more than compensates for the effect of extra 
opacity. For example with an O abundance of log eq — 8.55, 
the predicted equivalent width of the OH 312.80nm line with 
the ID MARCS model atmosphere, ^tuib = l.Okms^^ and no 
extra opacity is 5.60pm, i.e. weaker than the "observed" line 
strength. In other words, this would at face value imply that 
there is no missing UV-opacity but rather that the known opac- 
ity sources predict too much opacity in this particular wave- 
length region. However, the reader is urged to be cautious and 
not over-interpreting the very preliminary non-LTE results of 
Asplund & Garcia Perez (2001). It is likely that the inclusion 
of more levels in the OH model would more efficiently ther- 
malize the molecule, bringing the results closer to LTE. It is 
clear that a more careful non-LTE investigation of the OH line 
formation is necessary to resolve this issue. 

Non-equilibrium chemistry for OH molecule formation: 
Another crucial assumption for these types of studies is the use 
of instantaneous equilibrium chemistry for the molecule for- 
mation. In reality the time-scale for molecule formation is fi- 
nite and if this approaches the typical convective time-scale 
for up- and down-flows in the solar atmospheres the number 
density of molecules will in general be lower than expected 
from equilibrium chemistry. This has been shown by Asensio 
Ramos et al. (2003) for the case of the CO molecule in hy- 
drodynamical models of the solar chromosphere. As for non- 
LTE line formation for molecules, little work has been done 
in this potentially crucial area. The available evidence, how- 
ever, indicates that equilibrium chemistry is valid for OH, at 
least in the atmospheric layers corresponding to the forma- 
tion heights of the OH vibration-rotation and A-X lines in the 
Sun (Sanchez Almeida et al. 2001; Asensio Ramos et al. 2003; 
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Asensio Ramos & Trujillo Bueno 2003). In addition, these ef- 
fects only influence the present study in as far as the two types 
of lines are formed in different layers with a differing degree of 
appropriateness for the assumption of equilibrium chemistry: if 
both the vibration-rotation and A-X lines have the same typi- 
cal formation heights any non-equilibrium molecule formation 
will affect both similarly, leaving the estimate of the missing 
UV-opacity and the solar Be abundance intact. 

Non-LTE line formation for Be lines: Finally, possi- 
ble departures from LTE for the line formation of the Bell 
lines can of course influence the derived solar Be abundance. 
Fortunately, the two main non-LTE effects for these lines - 
over-ionization and over-excitation - essentially compensate 
each other for the solar case, coincidentaUy leaving the non- 
LTE Be abundance very similar to the LTE value: A(log eBo) < 
0.02 dex (Garcia Perez et al., in preparation). Thus, the finding 
of no significant Be depletion is safe to departures from LTE in 
the Be line formation. 

5. Conclusions 

In the present article a re-analysis of the work by Balachandran 
& Bell (1998) has been performed employing a state-of-the- 
art 3D hydrodynamical solar model atmosphere (Asplund et al. 
2000a) instead of traditional ID hydrostatic models. Several of 
the uncertainties attached to the analysis by Balachandran & 
Bell (1998) are thus removed and the important results placed 
on a more firm footing. In addition, several of the outstand- 
ing sources of uncertainties have here been discussed at some 
length and, with one notable exception, found to be in most 
likelihood unimportant for the basic findings. The conclusions 
by Balachandran & Bell (1998) of a significant missing UV- 
opacity around 313 nm and a photospheric Be abundance in 
very close agreement with the meteoritic value measured in 
Cl-chondrites are thereforeconfirmed by this study based on 
a more realistic model atmosphere. According to this study, 
the additional opacity amounts to about 50% in this partic- 
ular wavelength region, which is only slightly smaller than 
Balachandran & Bell (1998) estimated using ID model atmo- 
spheres (about 60%). 

Although this study does not shed further light on the na- 
ture of the opacity source causing the missing UV-opacity it is 
likely that it can be traced to photo-ionization of Fe I. Using 
data from the Iron Project, Bell et al. (2001) indeed found that 
such unaccounted for Fe I bound-free cross-sections are close 
to explaining the missing 60% continuous opacity estimated by 
Balachandran & Bell (1998). This will ring even more true with 
the sUghtly smaller multiplication factors derived here. 

The situation related to the depletion of the light ele- 
ments in the solar photosphere can thus be summarized as fol- 
lows: while Uthium has been destroyed by a factor of about 
140 (Lodders 2003; Asplund et al. 2003a), neither the pho- 
tospheric berylHum nor the boron abundances appear to have 
been significantly modified throughout the solar evolution to 
the present day. As described in the Introduction, no standard 
model of stellar evolution can account for this depletion pat- 
tern. Furthermore, is poses a severe chaUenge for non-standard 
models such as slow rotational-induced mixing below the con- 



vection zone. Indeed, it is a benchmark which any such models 
must successfully be tested against before they can be trusted 
for other types of stars. 
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